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ABSTRACT 

We report on measurements of the cosmological constant, A, and the redshift space 
distortion parameter /3 = Q^/b, based on an analysis of the QSO power spectrum 
parallel and perpendicular to the observer's line of sight, P s (k\\ , kj_), from the final 
catalogue of the 2dF QSO Redshift Survey. We derive a joint A — j3 constraint from 
the geometric and redshift-space distortions in the power spectrum. By combining 
this result with a second constraint based on mass clustering evolution, we break 
this degeneracy and obtain strong constraints on both parameters. Assuming a flat 
(f2 m + J7a = 1) cosmology and a A cosmology r(z) function to convert from redshift 
into comoving distance, we find best fit values of ft a = 0-711^17 and M z ~ 1A ) = 
0.45^0 ^. Assuming instead an EdS cosmology r(z) we find that the best fit model 

obtained, with fl\ = 0.64+°'^ and (3 q (z ~ 1.4) = 0.40ig'og, is consistent with the 
A r(z) results, and inconsistent with a £1a = flat cosmology at over 95 per cent 
confidence. 

Key words: cosmology: observations, large-scale structure of Universe, 
quasars: general, surveys - quasars 



1 INTRODUCTION 

Observations of high redshift Type la Supernovae (SNIa) 
have recently generated much excitement and interest in 
ACDM cosmological models, deriving strong constraints on 
the cosmological constant, A (e.g. Perlmutter et al. 1999; 
Riess et al. 1998). There are still uncertainties regarding the 
standard candle assumption that underlies this approach, 
however, and we don't yet have a good understanding of 
SNIa physical processes. It is therefore important to develop 
independent ways to constrain A in order to confirm this 
tantalizing result, as different methods suffer from different 
systematic uncertainties. 

Different astronomical datasets are sensitive to different 
combinations of the cosmological parameters. For example, 
the recent WMAP cosmic microwave background (CMB) 
observations provide a strong constraint on £l m h 2 (Spergel 
et al. 2003), whereas the shape of the matter power spec- 
trum, probed by either galaxies (Percival et al. 2001) or 



QSOs (Outram et al. 2003), is sensitive to Q m h. By com- 
bining several approaches we can break these degeneracies 
and hence derive much stronger constraints on the model 
that best describes our Universe (e.g. Efstathiou et al. 2002). 
Neither the CMB nor the large-scale structure observations, 
however, are directly sensitive to Oa, and can only indirectly 
infer its value when the two datasets are combined. 

Alcock & Paczyhski (1979) suggested that A might 
be measured directly from redshift-space distortions in the 
shape of large-scale structure, by making the simple assump- 
tion that clustering in real-space is on average spherically 
symmetric. Geometric distortions occur if the wrong cos- 
mology is assumed, due to the different dependence on cos- 
mology of the redshift-distance relation along and across the 
line of sight, and from the size of these geometric distortions, 
the true cosmology can be determined. 

One approach to implementing this test is to use the 
distribution of QSOs to probe large-scale structure. By com- 
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paring the clustering of QSOs along and across the line of 
sight and modelling the effects of peculiar velocities and bulk 
motions in redshift space (Ballinger et al. 1996, Matsubara 
& Suto 1996), geometric distortions can be detected. 

The ideal sample for this is the recently completed 2dF 
QSO Redshift Survey (2QZ). Data for 2QZ were obtained 
using the AAT 2dF facility. The completed survey comprises 
some 23000 bj < 20.85 QSOs in two 5 x 75 deg 2 declination 
strips, one at the South Galactic Pole and one in an equato- 
rial region in the North Galactic Cap spanning the redshift 
range 0.3 z 2.5. The 2QZ catalogue and spectra were re- 
leased in July 2003 (Croom et al. 2003) and can be obtained 
at http://www.2dfquasar.org. 

Outram et al. (2001) produced a preliminary anal- 
ysis of redshift-space distortions in the power spectrum, 
P s (fc||,kx), of the incomplete 10k 2QZ survey (Croom et 
al. 2001). To help develop the method, and to test estima- 
tors designed to discriminate the effects of infall (measured 
via the parameter /3 ~ Q,^ 1 /b), small-scale velocity disper- 
sion, and geometric distortions, a simulation of the 2QZ was 
used. Three light cone strips of the redshift survey were sim- 
ulated, using the Virgo Consortium's huge Hubble Volume 
N-body ACDM simulation (Frenk et al. 2000, Evrard et al. 
2002). The effects of infall and geometry on the clustering 
distribution are very similar, and this approach alone could 
therefore only provide a degenerate constraint on A and (3. 
Outram et al. found, however, that this degeneracy could be 
broken by jointly considering a second constraint based on 
mass clustering evolution. Using the Hubble Volume sim- 
ulation, Outram et al. predicted that this method should 
constrain /3 to approximately ±0.1, and Qa to ±0.25 us- 
ing the final 2QZ catalogue. A complementary approach, 
using the two-point correlation function, £(<r, 7r), to inves- 
tigate redshift-space distortions in QSO clustering can also 
be considered (Hoyle et al. 2002). Here we shall apply the 
analysis developed in Outram et al. to the final 2QZ sample, 
and hence derive this significant new constraint on A. 



2 POWER SPECTRUM ANALYSIS 

We reapply the analysis described in Outram et al. (2001) to 
the final 2QZ catalogue, containing 22652 QSOs (only those 
QSOs with quality 1 are used in this analysis; see Croom 
et al. (2003) for further details). For convenience we review 
the main details below. 

2.1 The 2QZ selection function 

We account for the various selection effects introduced into 
the 2QZ catalogue by generating a catalogue of random 
points that mimics the angular and radial selection func- 
tions of the QSOs but otherwise is unclustered. Details of 
the QSO selection function used can be found in Section 
2.1 of Outram et al. (2003). We are measuring QSO clus- 
tering as a function of comoving distance, and so need to 
assume a cosmology in order to convert from redshift into 
comoving distance. Geometric distortions in the shape of 
large-scale structure occur if this cosmology is wrong, due 
to the different dependence on cosmology of the redshift- 
distance relation along and across the line of sight. From 
the size of these geometric distortions, measured in the QSO 
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Figure 1. P s (fc||,kj_) determined from the 2QZ Cata- 
logue assuming an EdS r(z). Filled contours of constant 
\og(P(k)/h~ 3 Mpc 3 ) are shown as a function of fcy/ZiMpc -1 and 
kx//iMpc . Overlaid are the best fit model (solid contours) 
with (3 = 0.40 and Qa = 0.64, and for comparison, the best fit- 
ting flat Qa = model with fi = 0.27 (dashed contours). 



power spectrum, P S (k\\ , kx), the true cosmology can be de- 
termined. The choice of adopted cosmology shouldn't matter 
in this analysis; adopting an incorrect cosmology does not 
affect the validity of this result. To check for consistency, 
however, we consider two possibilities; an Einstein-de Sitter 
(fi m =1.0, f2A=0.0) cosmology r(z) (EdS hereafter) and an 
f2 m =0.3, Qa=0.7 cosmology r(z) (A hereafter). To limit in- 
completeness, we restrict our analysis to 0.3 < z < 2.2. This 
restricts our sample to 19549 QSOs with a mean redshift of 
z~ 1.4. 



2.2 Power spectrum estimation 

The power spectrum estimation is carried out as described 
in Outram, Hoyle & Shanks (2001). When calculating 
P s (fc||,kx), information about the line of sight must be re- 
tained. To achieve this the data are divided into subsamples 
that subtend a small solid angle on the sky, and the distant- 
observer approximation is then applied. The data from each 
declination strip are split into 8 regions each approximately 
5° x 10° during this analysis. Each region is embedded into 
a larger cubical volume, which is rotated such that the cen- 
tral line of sight lies along the same axis of the cube in each 
case. The density field is binned onto a 256 3 mesh, using 
nearest grid-point assignment. The power spectrum of each 
region is estimated using a Fast Fourier Transform (FFT), 
and the average of the resulting power spectra is taken. The 
results, binned logarithmically into k\\ and kx, are plotted 
in Figure 1 (assuming an EdS r(z)) and Figure 2 (assuming 
a A r(z)). 
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Figure 2. P s (feii,kx) determined from the 2QZ Catalogue as- 
suming a A r(z). Filled contours of constant log(P(fc)/ft~ 3 Mpc 3 ) 
are shown as a function of fcy/ZiMpc -1 and k±/hVlpc~ 1 . Over- 
laid is the best fit model with fi = 0.45 and SIa = 0.71. 



2.3 The window function 

The measured power spectrum is convolved with the power 
spectrum of the window function. The power spectrum of the 
window function is shown in figure 3. We compare the 2QZ 
window function to that of the 10k catalogue and the 25k 
mock catalogue determined in Outram et al. (2001). Due to 
the varying observational and spectroscopic incompleteness 
imprinted on the angular distribution of the 2QZ QSOs, the 
points lie slightly above those of the mock catalogue (where 
it was assumed that the angular distribution of the final 
catalogue would be uniform), but significantly lower than 
those of the 10k catalogue, due to the increased coverage. 
Outram et al. (2001; Section 2.3) adopted limits of fcy & 
|kj_| > 0.02 h Mpc -1 for the analysis of P s (fc||,k ± ) to avoid 
any problems due to the effect of the window function on 
large scales. Whilst slightly higher than that of the mock 
catalogue, the window function is still sharply peaked in k- 
space, and so we adopt the same limits in the case of the 
A r(z) analysis. Any window function effects should still be 
negligible. Assuming an EdS r(z), however, the comoving 
distances probed are smaller, and hence the window function 
affects smaller scales. The appropriate limits to adopt in this 
case are kn & |kj_| > 0.03 h Mpc~\ 



2.4 Error estimates 

A different error estimate is adopted in this paper from that 
of Outram et al. (2001). They adopted two methods to esti- 
mate the error in the power spectrum measurements; either 
via the dispersion between the power spectrum measure- 
ments from three declination strips of the Hubble Volume 
simulation, or the error estimate obtained using the method 
of Fcldman, Kaiser & Peacock (FKP; 1994). The two es- 
timates agree fairly well, but with large scatter. The FKP 
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Figure 3. The power spectrum of the window function of the 2QZ 
catalogue, assuming A r(z), shown as a function of k. The solid 
squares show the points which were used in the analysis, with 
wavenumbers with k\\ > 0.02, |k x | > 0.02, and k < 0.2. The dots 
show the power spectrum at smaller wavenumbers; excluded from 
this analysis. The window function is a steep power law, varying as 
~ k~ 3 . For comparison, the open circles show the power spectrum 
of the window function from the analysis of the 2QZ 25k mock 
catalogue, and the crosses show the power spectrum of the window 
function from the analysis of the 2QZ 10k catalogue (Outram et 
al. 2001). The points are significantly lower than those of the 10k 
catalouge, due to the increased coverage, but still slightly higher 
than those of the mock catalogue due to the varying completeness 
pattern that still affects the finished survey (Croom et al. 2003). 
The overlaid line is a fit to the window function from the usual 
power spectrum analysis, binning in shells of k (Outram et al. 
2003). 



error is, on average, slightly lower than the error estimated 
from the dispersion betwen the three mock strips. Neither 
approach is ideal; the error estimated from the dispersion 
has a high uncertainty due to the small number of mock 
catalogues, and relies on the power of the mock catalogues 
being similar to that of the 2QZ catalogue, whereas the FKP 
error estimate assumes a thin shell geometry, and can under- 
estimate the error on some scales, due, for example, to the 
onset of non-linearity. 

In order to apply the distant-observer approximation 
(see Section 2.2) the QSO power spectrum was estimated 
separately in sixteen individual 5° x 10° regions. For this 
paper we have chosen to estimate the errors self-consistently 
from the data by considering the dispersion in these sixteen 
power spectrum measurements. As these regions are taken 
from just two contiguous strips, they lie next to each other, 
and are not entirely independent. It is therefore possible that 
they slightly underestimate the error in the power estimate, 
especially on large scales. Whilst we believe this effect would 
be small on the scales of interest, and hence do not consider 
it further in this paper, to test this fully would require sev- 
eral independent simulations of the survey, which are not 
currently available due to the large volume probed. This er- 
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ror estimate agrees well, on average, with estimates obtained 
from the Hubble Volume simulation (see section 2.4 of Out- 
ram et al. 2001), but is slightly larger than the FKP error 
estimate. 



3 MODELLING REDSHIFT-SPACE 
DISTORTIONS 

The power spectrum model incorporating redshift distor- 
tions that we fit to the data is presented in Ballinger, Pea- 
cock and Heavens (1996), and briefly summarized below. 

The power spectrum analysis was carried out assum- 
ing either an fi m =1.0, fiA=0.0 cosmology or an f2 m =0.3, 
riA=0.7 cosmology. If the true cosmology differs from this 
then geometric distortions will be introduced into the clus- 
tering pattern, due to the different dependence on cosmology 
of the redshift-distance relation along and across the line of 
sight. Our distance calculations will be wrong by a factor 
f± perpendicular to the line of sight, and fu along the line 
of sight (as defined in Ballinger et al. 1996). Thus we can 
define a geometric flattening factor: 



(1) 



(Ballinger et al. 1996). The effect this has on the power 
spectrum is given by 



P ic(fc||,k x ) = 



■ft rue (&) 



f 3+n p 



(2) 



(Ballinger et al. 1996) where jx = k\\/k, and n is the spec- 
tral index of the power spectrum. By measuring the size 
of this geometric flattening, F, seen in the power spectrum 
P s (fe||, kx), we can therefore calculate the true value of SIa- 
Unfortunately the problem is complicated by the fact 
that redshift-space distortions are also caused by peculiar 
velocities. The main cause of redshift-space distortions on 
the large linear scales probed by the QSO power spectrum 
are coherent peculiar velocities due to the infall of galaxies 
into overdense regions. This anisotropy takes a very simple 
form in redshift-space, depending only on the density and 
bias parameters via the combination (3 « Q, 1 ^ 6 /b: 

21 2 



P S (fc|,,kx) = P R (k) [l + /3/J 2 ]' 



(3) 



(Kaiser 1987) where P s , and P R refer to the redshift-space 
and real-space power spectra respectively. 

Whilst we expect the parameters b, f2 m , and hence (3 
to vary as a function of redshift, we determine an average 
P s (fe|j,kx) over the redshift range 0.3 < z < 2.2 from the 
2QZ data, and hence only fit a single value of f3 q {z ~ 1.4) to 
the observed distortion. We consider the effect of allowing 
P q (z) to vary on the model, to test whether such variations 
could introduce systematic uncertainties, and whether the 
value of (3 obtained represents a true average value for (3 q 
over the range 0.3 < z < 2.2. Combining the redshift-space 
distortion power spectrum models generated using a varying 
f3 q (z), weighted using the QSO n(z), we can determine an 
average model P s (fc| h kx). We assume a simple bias model, 
designed to match the observations of Outram et al. (2003) , 
where the QSO clustering amplitude increases linearly with 
redshift by 30 per cent from z = 0.3 to z = 2.2 (assuming A). 
In this evolving bias model, the best-fitting value of (3 q (z) 



increases with redshift from /3 q (0.3) = 0.45, reaching a maxi- 
mum of /3,(0.75) = 0.51, before decreasing to (3 q {lA) = 0.45 
and (3 q {2.2) — 0.34. The value averaged over the QSO n(z) 
is /3 q = 0.44, in almost perfect agreement with the simple 
model with a single value of f3 q = 0.45, and the value of 51a 
obtained was unchanged, indicating that any systematics in- 
troduced by variations in /3 with redshift are much smaller 
than the errors. 

To truly constrain {3 as a function of redshift, using 
redshift-space distortions, we would have to measure the 
QSO power spectrum in redshift bins. Unfortunately, due 
to the smaller sample sizes and volume probed for each bin, 
the errors in such an analysis would be considerably larger 
than any expected change in (3 with the size of the cur- 
rent dataset, and so no useful constraint would be obtained. 
Therefore we have not included such an analysis in this pa- 
per. 

Anisotropy in the power-spectrum on small scales is 
dominated by the galaxy velocity dispersions in virialized 
clusters. This is modelled by introducing a damping term. 
The line of sight pairwise velocity (<r p )* distribution is mod- 
elled using a Lorentzian factor in redshift-space: 



D [kfMTp] 



1 + 5 (kfia p ) 



(4) 



Combining these effects leads to the final model: 



p s (W = P {k) 



fl +n F 

1 + ^ 2 (^ 1 - 1 )] 2 - D [ fc K] 



(5) 



(Ballinger et al. 1996) where a' p = <r p /f\\ ■ 



4 MASS CLUSTERING EVOLUTION 

Following Outram et al. (2001), we also consider a second 
constraint on A and (3, based on the evolution of mass clus- 
tering from the average redshift of the 2QZ to the present 
day. 

First we need an estimate of the z ~ mass cluster- 
ing amplitude. For this we could adopt local measurements 
of as (e.g. Eke et al 1996; Hoekstra et al 2002), however, 
to limit systematics, we choose to determine the local mass 
clustering amplitude using the reverse of the method ap- 
plied at z ~ 1.4 using 2QZ. Therefore we will use the recent 
observations of the 2dF Galaxy Redshift Survey (2dFGRS). 

Hawkins et al. (2003) used the 2dFGRS catalogue to 
determine the real-space galaxy correlation function at an 
effective redshift of z = 0.15, finding it to be well-described 
by a power-law (r/ro) _7r out to scales of r ~ 20fo~ 1 Mpc, 
where ro = 5.05/i -1 Mpc and j r — 1.67. Using the same 
sample of galaxies, Hawkins et al. produced a measurment 
of /3 g (z ~ 0) = 0.49 ± 0.09, using redshift-space distortions 
in the two-point galaxy correlation function. The value of 
f3 g was measured on scales 8 < s < 30ft _1 Mpc. 

For each cosmology (again we only consider flat cos- 
mologies) the value of the galaxy-mass bias can be found 



* In power spectra a p is implicitly divided by Ho and quoted in 
units h- 1 Mpc. (H = WOh km s" 1 Mpc" 1 ) 
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from P s , which in turn gives the amplitude of the mass cor- 
relation function from the measured galaxy correlation func- 
tion. The evolution in the mass correlation function from 
z = 0.15 to z = 1.4 can then be calculated for this cosmol- 
ogy using linear theory, and hence by comparing the z — 1.4 
mass correlation function to the 2QZ correlation function, an 
estimate of the QSO bias factor, and therefore f3 q (z ~ 1.4) 
can be obtained, as a function of cosmology. As in the pre- 
vious section, this method assumes that a single value of f3 q 
applies at the average redshift of the QSO survey, which, 
at least in the case of the evolving bias model discussed in 
Section 3, is a good approximation (we obtain a best fit of 
P q (z = 1.4) = 0.45 in the evolving bias model, assuming A, 
almost identical to the average value of q = 0.44). 



5 RESULTS 

5.1 Fitting the Redshift-space Distortions 

Here we fit the redshift-space distortions model discussed 
in Section 3 to the QSO power spectrum, P s (fey,kx), de- 
scribed in Section 2. There are several free parameters in 
the model (Equation 5); F,f3,a p , and P R (k), the underly- 
ing real-space power spectrum. To reduce the uncertainties 
on the parameters of interest, /3 and F, we make simple 
assumptions about the other parameters. There is an uncer- 
tainty in determining QSO redshifts from low S /N spectra of 
8z ~ 0.0038 (Croom et al. 2003). This introduces an appar- 
ent velocity dispersion of a p ~ 670 km s _1 . By adding this 
in quadrature to the intrinsic small-scale velocity dispersion 
in the QSO population, we therefore expect an observed ve- 
locity dispersion of a p ~ 670 — 840 km s _1 (assuming an 
intrinsic velocity dispersion of a v ~ — 500 km s _1 ). a p is 
not well constrained in this analysis, due to the large scales 
probed and so we have allowed a p to vary, within the range 
670 < <7 P < 840 km s _1 , choosing the value that maximised 
the likelihood of the fit. Whilst a p has relatively little effect 
on the power spectrum at large scales, any difference in the 
true value from our assumed range would lead to a small 
systematic shift in the resulting best fit values of 51a and 
p. This could be due, for example, to an underestimate of 
the QSO redshift determination error caused by the intrinsic 
variation in line centroids between QSOs. 

We reconstruct a real-space power spectrum self- 
consistently from the measured QSO redshift-space power 
spectrum (Outram et al. 2003). To remove noise we adopt 
the best-fitting model redshift-space power spectra given in 
Table 4 of Outram et al.. The real-space power spectrum 
is then obtained by inverting the above redshift distortion 
equations in each cosmology for an average value of /j. Fi- 
nally, we assume a flat cosmology, consistent with the recent 
WMAP CMB results (Spergel et al. 2003). We choose to fit 
the variable Qa, and fix fi m = 1 — S7a- 

We fit the model to the 2QZ data by performing a maxi- 
mum likelihood analysis, allowing the parameters /3, and £Ia 
to vary freely and determining the likelihood of each model 
by calculating the \ 2 value of each fit. Only those wavenum- 
bers with fey > 0.02, |kj_| > 0.02, and k < 0.2 (assuming A), 
or fey > 0.03, |kj_| > 0.03, and k < 0.3 (assuming EdS) are 
used in the fit. The former constraints are applied to remove 
the effects of the window function, and the latter because the 
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Figure 4. Likelihood contours are plotted in the S7a - /3 plane 
for x 2 values corresponding to a one-parameter confidence of 68 
per cent, and two-parameter confidence of 68 and 95 per cent 
(dashed contours), calibrated using Monte Carlo simulations, for 
fits to P s (fc||,kx) determined from the 2QZ catalogue assum- 
ing an EdS r(z). Overlaid are the best-fit (dot-dash) and 1-a 
(dot) values of /3 determined using the mass clustering evolution 
method. Significance contours given by joint consideration of the 
two constraints are also plotted for a one-parameter confidence of 
68 per cent, and two-parameter confidence of 68 and 95 per cent 
(solid contours). The best fit model obtained (marked with a +) 
has = 0.40 and Q A = 0.64. 

FFT is unreliable at smaller scales. This constraint also pre- 
vents excessive non-linearity, where the model breaks down. 

Figures 4 and 5 show likelihood contours in the 51a 
- j3 plane, assuming an EdS and A cosmology r(z) re- 
spectively. Nominally the best fit values obtained from the 
redshift-space distortions analysis are f3 = 0.301q^9, and 
Qa = 0.76iQ'3g, w ilh x 2 — 61.5 over 57 degrees of freedom 
assuming EdS, or /3 = OA3±%% 9 , and Q A = 0.73±g;|g , with 
X 2 = 54.1 over 61 degrees of freedom assuming A. In both 
cases, the absolute \ 2 values indicate that the best-fitting 
models provide an adequate fit to the power spectrum. How- 
ever, when comparing the absolute x 2 values from different 
power spectrum realisations, we caution that there is consid- 
erable (~ 10 per cent) noise in the \ 2 values obtained, which 
manifests itself even with the same cosmology and selection 
function, when only subtle changes to the FFT parameters 
or binning are made. As Figures 4 and 5 show, there is a large 
degeneracy between the fitted values of the two parameters, 
due to the similarity in the shape of the redshift-space and 
geometric distortions (Ballinger et al. 1996). 

5.2 Testing the likelihood contours 

To test and calibrate the likelihood contours and hence the 
parameter uncertainties derived from the \ model fitting 
to the 2QZ P s (fc||,kjJ, Monte Carlo simulations were per- 
formed. 1000 realizations of the A r(z) 2QZ power spec- 
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Figure 5. Likelihood contours are plotted in the Qa - /3 plane 
for x 2 values corresponding to a one-parameter confidence of 68 
per cent, and two-parameter confidence of 68 and 95 per cent 
(dashed contours), calibrated using Monte Carlo simulations, for 
fits to P s (fc||,kx) determined from the 2QZ catalogue assuming 
a A r(z). Overlaid are the best-fit (dot-dash) and 1-a (dot) val- 
ues of (3 determined using the mass clustering evolution method. 
Significance contours given by joint consideration of the two con- 
straints are also plotted for a one-parameter confidence of 68 per 
cent, and two-parameter confidence of 68 and 95 per cent (solid 
contours). The best fit model obtained (marked with a +) has 
P = 0.45 and Q A = 0.71. 
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Figure 6. A comparison between the two-parameter confidence 
of 68 and 95 per cent likelihood contours (dashed), derived us- 
ing the x 2 statistic from a fit to P s (fc|| , kj_) measured from the 
2QZ catalogue assuming a A r(z), and solid contours derived from 
fOOO Monte Carlo simulations of P s (fc|| , kx) drawn from the best 
fit model (marked with a +), assuming the fractional errors on 
the power spectrum calculated from the dispersion in the mea- 
surements from 5° X f0° regions. The points show the best fit 
values for Hyy and obtained from each realization of the power 
spectrum. The straight solid and dashed lines show the median 
and f6/84 percentile values for Qa and j3 from the Monte Carlo 
simulations, which can be compared to the inner dashed contour 
showing the one-parameter 68 per cent confidence likelihood con- 
tour. 



trum were drawn from the best fitting model (/3 = 0.45, 
Qa = 0.71), assuming a fractional uncertainty on each data 
point equal to that estimated from the true power spectum 
(see Section 2.4). The redshift-space distortions model was 
fitted to each realization of the QSO power spectrum, as dis- 
cussed above. The best fit values for Qa and f3 obtained from 
each realization of the power spectrum are plotted in Fig. 6, 
together with the median and 16/84 percentile values of 
each parameter. The best fit values obtained from the Monte 
Carlo redshift-space distortions analysis are (3 — 0. 451^32, 
and Qa = 0.70_"q'4q, in good agreement with the input model 
to the simulations (f3 = 0.45, Qa = 0.71), however, the size 
of the uncertainties estimated via the Monte Carlo analysis 
are slightly larger than those estimated via the \ 2 likeli- 
hood analysis. The \ 2 contours containing 68 and 95 per 
cent of the simulated measured parameters were also calcu- 
lated from the simulations, and are compared to the equiva- 
lent % 2 likelihood contours in Fig. 6. Again, the Monte Carlo 
contours appear fairly close, but slightly larger than the like- 
lihood contours. The Monte Carlo simulations indicate that 
the likelihood contours underestimate the values corre- 
sponding to a given confidence level by approximately 20 per 
cent, and so for this analysis we have corrected the contour 
levels by this factor. 



5.3 The Mass Clustering Evolution Constraint 

We now combine the degenerate redshift-space distortion 
constraint obtained from the likelihood analysis above with 
the constraint obtained from consideration of the evolu- 
tion of mass clustering, described in Section 4. The value 
of /3 q (z ~ 1.4) and the one-sigma errors from this constraint 
are plotted on figures 4 and 5. Although totally degenerate 
with the value of Qa, this provides a different, and almost 
orthogonal constraint to that from the redshift-space dis- 
tortions. This was derived using significantly smaller scales 
than the P s (k\\ , kx) analysis, and so we can treat the results 
as independent and hence combine the likelihoods, yielding 
a much stronger fit. 

Assuming an EdS r(z), the joint best fit values obtained 
are Qa = 0.641°;^ and /3 q (z ~ 1.4) = 0.40_"g;°g. The best 
fitting flat Qa = model with (3 = 0.27 (shown in Figure 1) 
can be excluded at over 95 per cent confidence. The joint 
best fit model to the A r(z) 2QZ data has Q A = 0.7lt%°r 
and /3 q (z ~ 1.4) = 0.45j^J n- Likelihood contours in the Oa 
- (3 plane for the two fits are plotted for a one-parameter 
confidence of 68 per cent, and two-parameter confidence of 
68 and 95 per cent in figures 4 and 5. 

By combining this constraint with that derived from 
the redshift-space distortions, we are comparing values of 
(5 that were measured using different estimators and on 
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different scales; r ~ 20h 1 Mpc for the mass clustering 
evolution method described in this section, compared to 
r ~ 100/i _1 Mpc for the redshift -space distortions in the 
power spectrum. Hence we are implicitly assuming that bias 
is scale independent on these scales. Whilst this is likely to 
be true (e.g. Verde et al. 2002), a non- linear bias on these 
scales would introduce a systematic error in our results. To 
reduce any possible systematic effects further the value of (5 g 
could be measured using the same method (Outram, Hoyle 
& Shanks 2001), however this analysis has yet to be done us- 
ing the large 2dF galaxy sample. Although at a much lower 
redshift, there will also still be a small cosmological depen- 
dence in the determined value of (3 g which should be taken 
into account. 

Instead of using the correlation function amplitude to 
trace the evolution in clustering, the amplitude of the spher- 
ically averaged QSO power spectrum (Outram et al. 2003) 
could be compared to that of local galaxies (Percival et al. 
2001). The different effects on the respective power spec- 
tra of small-scale velocity dispersions, redshift uncertainties 
and, most importantly, the two survey window functions 
have to be taken into account in the comparison. The pub- 
lished 2dFGRS power spectrum was measured assuming the 
A cosmology r(z), so comparing its amplitude to that of the 
A r(z) QSO power spectrum, and also assuming, as before, 
/3 g (z ~ 0) = 0.49 ± 0.09 we would predict j3 q = 0.43 ± 0.09, 
in good agreement with the prediction of fi q = 0.45 ± 0.09 
obtained using the correlation function analysis. 

Alternatively, we could choose to constrain the local 

mass clustering amplitude using the relation as = (0.52 + 

/ - 0.04)n-°' 52+0 ' 13n '™ (Eke et al. 1996), determined from 

the evolution of rich galaxy clusters, instead of using the 

2dFGRS observations. With this approach we again obtain 

very similar results, with Q,a = 0.75lto'i7 an d Pq(z ~ 1.4) = 
n41 +o.09 



6 COMPARISON WITH THE HUBBLE 
VOLUME SIMULATION 

Outram et al. (2001) produced an analysis of redshift-spacc 
distortions in the power spectrum, P s (fc||,kx), of a sim- 
ulation of the 2QZ, created using the Virgo Consortium's 
huge Hubble Volume N-body ACDM simulation (Frenk et al. 
2000, Evrard et al. 2002). The parameters of the simulation 
are O b =0.04, fi C DM=0.26, fi A =0.7, H o =70 km s^Mpc" 1 
and the normalisation, <T8, is 0.9. The simulation can be 
used to produce three light-cone 75x5 degree declination 
strips extending to z ~4. A simple biasing prescription and 
sparse sampling is used to give the mass particles a similar 
clustering pattern and selection function as the final 2QZ 
survey. 

As we have refined the method of analysis for this pa- 
per, for a full comparison with the results presented here 
we repeat the analysis of the Hubble Volume simulation, us- 
ing exactly the same method as for the final 2QZ catalogue. 
Two 5 x 75deg 2 declination strips, each containing 12500 
mock QSOs, were used in the analysis; see Outram et al. 
for further details. To match the 2QZ analysis, the mock 
QSO redshifts were degraded to the noise level from 2QZ 
(Sz ~ 0.0038, corresponding to an apparent velocity disper- 
sion of <7 P ~ 670km s _1 ). As before, a v was allowed to vary, 




0.2 - 



0.0 . , " :'l : ' l 

0.0 0.2 0.4 0.6 0.8 1.0 
/J(z~1.4) 

Figure 7. Contours are plotted in the Qa - P plane for a one- 
parameter confidence of 68 per cent, and two-parameter con- 
fidence of 68 and 95 per cent (dashed contours) for fits to 
P s (ku,\i±) determined from the 25k Hubble Volume mock QSO 
catalogue assuming a A r(z). Overlaid are the best-fit (dot-dash) 
and 1-<t (dot) values of /3 determined using the mass clustering 
evolution method. Significance contours given by joint consider- 
ation of the two constraints are also plotted for a one-parameter 
confidence of 68 per cent, and two-parameter confidence of 68 and 
95 per cent (solid contours) . The best fit model obtained (marked 
with a +) has = 0.36 and Q A = 0.72. 



within the range 670 < a p < 840 km s _1 , choosing the value 
that maximised the likelihood of the fit. 

The results are shown in Fig. 7. The joint best fit model 
to the A r(z) mock QSO data has f2 A = 0.72±g"_5| and 
P q {z ~ 1.4) = 0. 36i^o7- These results are consistent both 
with the input values for the simulation, f3 q = 0.35 and 
Qa = 0.7, and the best fitting values obtained by Outram 
et al. (2001). The uncertainty in the estimates of the two 
parameters, however, are considerably smaller than those 
quoted by Outram et al.. This is because they incorrectly 
calculated the one-parameter uncertainties for /3 and £Ia 
from the two-parameter % 2 contours, significantly overes- 
timating the true uncertainty. The two-parameter 68 and 
95 per cent confidence levels, however, are very similar to 
those derived by Outram et al.. Using the Hubble Volume 
simulation, and the uncertainty estimates derived from the 
one-parameter 68 per cent likelihood contour, we would in- 
stead predict that this method should constrain f3 to ap- 
proximately ±0.07, and S7a to ±0.1 using the final 2QZ cat- 
alogue. 

In the 2QZ catalogue there remains a small level of in- 
completeness that was not considered in the simulation, and 
therefore there are approximately 2000 fewer QSOs. The un- 
certainties inherent in constructing a completeness mask to 
correct for this in the 2QZ will also add a small amount of 
noise that was not considered in the simulation, although 
Outram et al. (2003) demonstrated that any potential sys- 
tematics due to this are significantly smaller than the statis- 
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tical errors in the power spectrum determination. Hence, the 
uncertainties quoted above are slightly optimistic. However, 
the robust way that the Hubble Volume analysis returned 
almost exactly the true input parameters, coupled with re- 
markable similarity between both the results and the uncer- 
tainties from the simulation, shown in Fig. 7, and the data, 
shown in Fig. 5, gives us great confidence in the results we 
present here from the 2QZ catalogue. 



7 DISCUSSION AND CONCLUSIONS 

In this paper we have reported on measurements of the cos- 
mological constant, A, and the redshift space distortion pa- 
rameter (3 — Q^/b, based on an analysis of the QSO power 
spectrum parallel and perpendicular to the observer's line 
of sight, P s (fe||,kx), from the final catalogue of the 2dF 
QSO Redshift Survey. We have derived a joint A — f3 con- 
straint from the geometric and redshift-space distortions in 
the power spectrum. By combining this result with a second 
constraint based on mass clustering evolution, obtained by 
comparing the clustering amplitude of 2QZ QSOs at z ~ 1.4 
with that of 2dFGRS galaxies at z ~ (Hawkins et al. 2003) , 
we have broken this degeneracy, obtaining strong constraints 
on both parameters. 

Assuming a flat (fi m + £7 a = 1) cosmology and a A 
cosmology r(z) we find best fit values of Q,a = 0.711q'j7 
and j3 q (z ~ 1.4) = 0.45i°;°j. Assuming instead an EdS cos- 
mology r(z) we find that the best fit model obtained, with 
Qa = 0.64^°;^ and /3,{z ~ 1.4) = 0.40±g;g| , is consistent 
with the A r(z) results, and again strongly favours a A dom- 
inated cosmology. Indeed, the EdS cosmology r(z) result is 
inconsistent with a flat Qa = cosmology at over 95 per 
cent confidence. 

We determine the 2QZ power spectrum in a single 
adopted cosmology. Following Alcock & Paczyhski (1979), 
Ballinger et al. (1996) and Outram et al. (2001), we fit the 
size of the geometric distortions in the power spectrum that 
arise if this cosmology is incorrect, hence deriving a con- 
straint on the true value of £Ia- An alternative approach 
would be to measure the 2QZ power spectrum in each of 
the cosmologies considered in this paper (with < Qa < 1), 
rather than just an adopted cosmology, then measure the 
goodness-of-fit of the redshift-space distortions model (with- 
out geometric distortions) in that cosmology, determining 
which cosmology provides the best overall fit. Obviously, this 
would be expected to occur when there were no geometric 
distortions in the power spectrum. A comparison of the \ 2 
values obtained for the best fits to the A and EdS cosmology 
2QZ power spectra suggests that this method might provide 
an even tighter constraint on Qa- However, there is consid- 
erable (~ 10 per cent) noise in the absolute \ 2 value of the 
best fit to any power spectrum realisation (which manifests 
itself even when only subtle changes to the FFT parame- 
ters or binning are made, let alone changes in the adopted 
cosmology) . Whilst this has little effect on the best fitting j3 
and Q,a parameters in a given realisation (much smaller than 
the quoted uncertainties), it does inhibit the comparison 
of x 2 fits from different realisations. This approach would 
no doubt also confirm that A dominated cosmologies are 
strongly favoured by the 2QZ power spectrum, however, a 
much noisier picture would therefore emerge. 



Outram et al. (2003) presented an analysis of the 
spherically-averaged 2QZ power spectrum. They found very 
good agreement between the QSO power spectrum and the 
Hubble Volume ACDM simulation over the whole range of 
scales considered. From the shape of the QSO power spec- 
trum Outram et al. derived a low value of Q m h — 0.19±0.05 
(where h is the Hubble parameter) that is hard to reconcile 
with a standard (fi m = 1) CDM model. Their results are 
entirely consistent with the analysis described in this paper, 
of redshift-space distortions in the same QSO power spec- 
trum, and taken together they provide strong evidence in 
favour of the ACDM scenario. 
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